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ABSTRACT 


Steady magnetic field measurements of magnitude 30 to 100 gamma on the 

lunar surface impose problems of interpretation when coupled with the non- _ 

detectability of a lunar field at 0. 4 lunar radius altitude and the limb induced 
perturbations of the solar wind reported by Mihalov et al. at the Explorer orbit. 
The lunar time varying magnetic field clearly indicates the presence of eddy 
currents in the lunar interior and allows calculation of an electrical conductivity 
profile. The problem is complicated by the day-night asymmetry of the moon's 

y* 

electromagnetic environment, the possible presence of the TM mode and the 
variable wave directions of the driving function. The electrical conductivity is 
calculated to be low near the surface, rising to a peak of 6 x 10 -3 ohm -1 meter -1 
at 250 km, dropping steeply inwards to a value of about 10 -5 ohm -1 meter -1 , and 

then rising toward the interior. A transition at 250 km depth from a high con- 

/ 

ductivity to a low conductivity material is inferred, suggesting an olivine-like 
core at approximately 800°C, although other models are possible. 


Introduction 

The measurement of lunar magnetic fields has long been of interest because 
of the promise that these measurements would provide information about the 
lunar interior: a significant dipolar field would imply a dynamo action in the 
lunar interior similar to the earth's dynamo, while the presence of induced fields 
would imply a value for the electrical conductivity in the lunar interior. 
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The first magnetic experiments were made by the USSR on the Luna space- 
craft. These set an upper limit on any lunar magnetic field and also indicated a 
magnetospheric region surrounding the moon. More definitive measurements of 
the lunar environment came with the launching of the Explorer 35 lunar orbiting 
satellite in 1967. Magnetometers on this satellite detected no lunar bow shock 

and no evidence at periselene of a lunar field (Colburn et al. , 1967, Ness et al. , 

20 

1967). The upper limit for a lunar centered dipole moment was set at 10 gauss 

q 

cm corresponding to a maximum dipolar field component of 4 gamma (1 gamma 
= 10“^ gauss) on the lunar surface (Behannon, 1968). Larger surface fields 
were not ruled out if they were of quadrupole or higher order. Definitive meas- 
urements of lunar fields commenced with the deployment of the Lunar Surface 
Magnetometer on Apollo 12 and the Lunar Portable Magnetometer on Apollo 14. 

In this paper we discuss first the measurement of a steady magnetic field 
at three points on the lunar surface. We consider next the predictions of other 
lunar magnetic regions inferred from satellite observations. Finally we con- 
sider the calculations of the electrical conductivity profile of the lunar interior 
determined by the moon’s electromagnetic response, and its implications for a 
thermal model. 

The Instruments 

Explorer 35 was placed into lunar orbit on July 19, 1967 and is still in 
operation. The periselene is at an altitude of 0.4 (R m = 1740 km) with the 
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aposelene altitude 4.4 R . The orbital plane is tilted 11° out of the ecliptic 
plane. Among the experiments are a plasma probe and two magnetometers. 

The Ames magnetometer that supplied the data reported here, is a three-axis 

✓ 

vector fluxgate with 0. 4 gamma (1 gamma _=A0 - Agauss) resolution- employing 

spin demodulation and filtering to avoid aliased data (Colburn et al. , 1967, 

Mihalov et al. , 1968). In this system the high bandwidth signals from the two 
sensors in the spin plane of the spinning spacecraft are multiplied by sines and 
cosines of the spin frequency and mixed appropriately to furnish a high bandwidth 
set of the three components of the vector field in an inertial (nonspinning) frame 
(Sonett, 1965). These are filtered to preserve the Nyquist criterion for alias free 
data before being sampled at the sampling rate (1 vector per 6. 14 seconds) and 
telemetered to earth (Sonett, 1968). The maintenance of spectral purity is vital 
to the use of power spectral density estimates in determining an electrical con- 
ductivity profile. 

The Apollo 12 Lunar Surface Magnetometer (LSM) is a triaxial fluxgate 
empiaced on the lunar surface at the Apollo 12 site at the eastern edge of Oceanus 
Procellarum, coordinates 3.0° S, 23.4° W. The instrument telemeters to earth 
a vector measurement of the field at a rate of 3. 3 samples per second. The site 
survey mode that was commanded three days after emplacement rotated the sen- 
sors to each of the three coordinate directions in turn. The three sensors are 
mounted at the ends of three mutually orthogonal 100 cm booms, permitting a 
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measurement of the local field gradient in the horizontal plane (Dyal et al. , 

1970b). Since no measurable gradient was found, a lower limit was placed on 
the distance to the field source, the distance also depending on the source con- 
figuration. A similar Lunar Surface Magnetometer was emplaced at the Apollo 15 
site, with a third one scheduled for Apollo 16. 

At the Apollo 14 site the field was measured by the Lunar Portable Magne- 
tometer (LPM). This instrument was set up at two locations by the astronauts 
who then sent back the readings by voice channel. A 20-second time constant 
filter averaged out the high frequency fluctuations to measure the steady state 
field component at the site. A similar Lunar Portable Magnetometer is scheduled 
to be operated on Apollo 16. 

Another magnetic measurement of the moon has become possible by the 
magnetometer aboard the subsatellite launched into lunar orbit by the Apollo 15 
command module (P. J. Coleman, Jr. , principal investigator). It is expected 
this instrument and the similar one on Apollo 16 will be able to map locally 
magnetized areas on the lunar surface when on the lunar dark side or in the 
earth's magnetic tail. The low orbital altitude, typically 100 km or less, will 
allow detection of surface feature at detail unobtainable at the orbit of Explorer 35. 
The Explorer measurement is very important in providing a measurement of the 
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background field surrounding the moon and is essentially unperturbed by the 
moon’s presence except when Explorer 35 is almost directly in the lunar shadow. 

— - - The~Surf ace Field- — — — 

Data from the Apollo 12 LSM showed a steady field component at the site of 
38 gamma and it was known at the outset that this was not due to a dipole centered 
in the moon, for a centered dipole of such strength would have been observed at 
the orbit of Explorer 35. It likewise could not have been due to a dipole source 
closer than 200 meters to the instrument because of the null result of the gradient 
measurement. The possibility of a localized surface field highly variable with 
location suggested magnetization of the lunar crustal regions and prompted the 
design and approval of the LPM that made its measurements only slightly more 
than one year after inception. 

Table 1 shows the vector field measurements reported for the two Apollo 14 
sites and the Apollo 12 site (Dyal et al. , 1971, Dyal et al. , 1970a). The vectors 

are also displayed in Fig. 1. The Apollo 12 and 14-C' (Cone Crater) magnitudes 

/ 

are approximately 40 gamma while the other Apollo 14 site has a magnitude factor 
of 2. 5 greater. The vectors point down and southerly; the maximum angle be- 
tween any two is 84°. The lower limit on the field gradient at Apollo 14 is de- 
fined as the average gradient between the two sites or 54 ±7 gamma/km, while 
the upper limit for Apollo 12, 133 gamma/km is determined by the site survey 


measurement, over a distance span on the order of 1 meter. 
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Table 1 


Magnetic Field Measurements at Lunar Sites (Dyal et al. , 1971) 


Site 

Apollo 12 

Apollo 14-A 

14-C' 

Coordinates 




S 

3.0° 

3.7° 

3.7° 

W 

23.4° 

17.5° 

17. 5 6 

Location 

LSM site 

170 m from LM 

Cone Crater Rim 

Field 

38 ±37 

103 ±57 

43 ±67 

Magnitude 




Components Up 

-24.4 ±2. 07 

-93 ±47 

-15 ±47 

East 

+13.0 ±1.87 

+38 ±57 

-36 ±57 

North 

-25. 6 ±0. 87 

-24 ±57 

-19 ±87 

Gradient 

<133 7/km 

>54 ±7 7/km 

• 


Models for the steady fields have been constructed but none have been com- 
pletely satisfactory. Since it is known from Explorer 35 measurements that 
the surface field caused by any centered dipole must be less than 4 gamma, more 
localized field sources are hypothesized, which would presumably rule out the 
mechanism of a planetary dynamo. Alternatively, the measurements could 
represent the bulk magnetization of surface rocks. The remanent field in a 
rock slab which has been uniformly magnetized has a highly nonuniform field, 
the magnitude being greatest near the edges and direction altering radically with 
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position near the edges. The magnitude becomes relatively smaller as one 
approaches the broad flat region of the magnetized slab. The edge delineated 
by Cone Crater, as pointed out by Dyal, or perhaps of a larger Mare region 
hidden by surface features (Runcorn and Quaide, private communication) would 
explain the large difference in the two Apollo 14 measurements. In view of 
these differences and the proposed models, the significance of all measurements 
lying in a single quadrant is not clear. More data would be required to establish 
the direction of the original magnetizing field. 

Remanent magnetism in rock areas consistent with magnetic measurements 
on Apollo rock can explain the steady field measurements (Runcorn et al. , 1970). 
This implies, however, that large rock volumes cooled through the Curie point 
under the influence of a field of more than 1000 gamma. This creates difficulty 
in lunar thermal history models, as will be discussed later. 

Magnetic Map Inferred from Explorer 35 

Mihalov et al. (1971) have found perturbations in Explorer 35 magnetometer 
measurements that appear to map the lunar surface. To describe the method, 
it is necessary to review the magnetic environment at the Explorer 35 orbit. 


The interplanetary magnetic field has a magnitude on the order of 5 gamma 
but can become several times greater under disturbed conditions. The preferred 



direction is outward or inward along an Archimedes spiral determined by the 
solar wind and the sun's rotation (Parker, 1963) although the field can assume 
any direction from time to time. The distinction between outward and inward 
can be made statistically, and for the period under question generally encom- 
passes one outward and one inward sector (occasionally more) every solar 
rotation period of 27 days (Wilcox and Colburn, 1970). The solar wind is not 
deflected significantly by the moon, as evidenced by the lack of a measurable 
bow wave, and hence the particles are believed to strike the lunar surface and 
be neutralized there, while directly behind the moon is a cavity in which solar 

J 

wind is essentially absent and the magnetic field is slightly enhanced because of 
the diamagnetism of the solar wind (Colburn et al. , 1971). At the border of the 
cavity a characteristic dip in magnetic field magnitude is observed. Beyond 
the Mach angle defining the rarefaction wave, however, a peak in the magnetic 
field is occasionally observed. These are believed to be due to a lunar solar 
wind interaction occurring at the limb and propagating outward along the Mach 
angle. The mechanism is postulated to be the deflection of solar wind by local 
magneticfieid. No magnetic effect of the moon is seen over the large majority 
of the Explorer 35 orbit lying forward of this Mach angle; the magnetic field is 
as if measured from a spacecraft very far from the moon. Because the solar 
wind travels at a speed faster than known magnetosonic and Alfvenic propagation 
speeds, the information signaling the presence of the moon cannot travel upstream 
to these regions. From time to time the magnetic field peak exterior to the 
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rarefaction wave becomes unusually large. There were 100 large peaks repre- 
senting 7.4% of the possible times of occurrence. Mihalov has catalogued these 
occurrences and related them to regions on the moon at the limb by associating 
each occurrence with the point on the lunar surface closest to the solar wind 
velocity vector drawn through the spacecraft location at the time of the observa- 
tion (Mihalov et al. , 1971, Sonett and Mihalov, 1971). The mapping has been 
shown to be statistically significant. In Table 2 are shown the areas outlined by 
Mihalov as being connected with the exterior peaks observed by Explorer 35. 

The table shows that the highlands are favored, and also that the far side is 
favored. While mechanisms other than local magnetic fields could be postulated 
as causes of the field increases seen by Explorer 35, Mihalov et al. find cogent 
reasons for believing that the local magnetic field is the most likely cause. On 
the theoretical side, Barnes et al. (1971) have investigated the possible frequency 
of magnetic structure on the surface . They find that regions of field such as 
indicated by the Explorer 35 perturbations must have a length scale of at least 
10 km and a compressed field strength of more than 10 gammas, and that thousands 
of these could be present on the lunar surface. 

Electrical Conductivity Profile Calculations 

The lunar surface magnetic field measurement was early thought of as a 
means of probing the interior of the moon (Sonett, 1966). If we consider the 
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. Table 2 

Regions on the lunar surface where the concentration of source locations was 
tested (Mihalov et al. , 1971). F is the fraction of observations showing the 
anomaly and P is calculated probability that there, is no concentration at the 
location. The observations are from the orbiting Explorer 35 satellite. Mihalov 
et al. conclude that local magnetization is the favored explanation for the 
phonomena. 


Seleno- 

centric 

latitude 

Seleno- 

centric 

longitude 

Area 

(X 10 5 km 2 ) 

F 

P 

Nearby future 

5°-20°S 

135°E-165°W 

10 

0.5 

«10 -5 

Mare X (Gagarin) 

6°N 

88° W 

2 

0.4 

<10“ 5 

"Montes d'Alembert" 

0°-20°N 

60°-115 o E 

10 

0.3 

3xl0 -3 

Mare Marginis 

5°S 

138°W 

0.7 

0.3 

6xl0 -3 

Crater 244 (Vavilov) 

2°N 

35°W 

1 

0.2 

<10" 5 

Eneke-Kepler 

5°N 

0° 

0.9 

0 . 1 

5xl0~ 4 

Pallas 

0° 

25°E 

0.8 

0.1 

4xl0 -5 

Delambre 


moon as a sphere whose electrical conductivity is a function of radius only, we 
find that there are two modes in which it can respond to the time varying elec- 
trical and magnetic field associated with the solar wind: the transverse magnetic 
(TM) and the transverse electric (TE) (Schubert and Schwartz, 1969). For the 
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TM response the driving function is an electric field due to the solar wind. 
Because of the high electrical conductivity of the plasma constituting the solar 
wind, the electric field must be considered to be zero in a frame at rest in the 
solar wind plasma. A simple transformation to a frame at rest with respect to 
the moon shows an electric field of E = v x B where B is the interplanetary 
magnetic field and v is the velocity of the moon with respect to the solar wind. 
This response should extend all the way down to zero frequency. The current 
system of the response must be such that the current flows through the moon 
and its crust and closes through the highly conducting solar wind. If this re- 
sponse were to be significant, it should cause a bow shock; a phenomenon not 
observed by Explorer 35. It is generally concluded that the relatively high 
resistance of the moon's crustal layers effectively cuts off the TM mode (Sonett 
and Colburn, 1967). The high resistance is deduced from the probable elec- 
trical conductivity of lunar material at the mean surface temperature, which 
must be assumed by material more than a few meters deep considering the low 
thermal conductivity of rock materials. An exhaustive statistical analysis of 
magnetometer measurements may yet show a modicum of TM response, but 
it must remain a secondary effect. 

We turn now to the TE mode, which is excited by B, the rate of change of 
the interplanetary magnetic field. The response of a planetary body to this 
mode is the establishment of eddy currents within the body that do not need to 
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flow through the relatively nonconducting crust. This mode does not occur in 
the steady state but is observable at frequencies commensurate with lunar 
magnetic field measurements. Preliminary findings on the TE response of the 
moon were reported by Sonett et al. (1971 b, c). 

Figure 2 shows a typical plot of power spectral density versus frequency 
for the horizontal components of the lunar surface field (Apollo 12) and for the 
same time period the corresponding components of the free stream solar wind 
field near the moon. The time period was 2 hours and the Apollo 12 magnetometer 
location 33° from the moon-sun line. The moon was outside the earth's bow shock 
in the solar wind. The power in the interplanetary field drops off with frequency 
approximately as f -2 , a typical frequency dependence for interplanetary spectra. 
The power in the surface components is larger, particularly at the higher 
frequencies. 

Autocorrelation techniques instead of cross correlation techniques are used 
because of the complication of the variable Doppler shift between the locations 
of the two measurements. The analysis is based on the assumption that the free 
stream field is the sole driving function for the lunar response. This is made 
credible by the examination of the amplification of particular sine wave cycles 
or transients such as were published in Sonett et al. (1971c). 
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Data from the first lunation were examined to find periods in which both 
magnetometers were reading simultaneously with no data gaps. Seven 2-hour 

l 

and seven 1-hour segments were used, with the Explorer 35 data transformed 
into LSM coordinates. The Apollo 12 data were filtered and decimated, to 
approximate the same filter and bit rate as the Explorer 35 data. Autocorrela- 
tion and power spectral density estimates were made for several frequencies. 

For the 2-hour swaths 50 frequencies were used, linearly spaced from 0. 83 
millihz to 41.5 millihz. For the 1-hour swaths the number of frequencies was 
25 (1.66 to 41. 5 MHz). The amplification A was defined as the square root of 
the ratio of power density at the lunar surface to power density in the free stream. 
The subscript x, y, or z is added to signify respectively the vertical, east, and 
north components of the vectors as viewed from Apollo 12. 


Values of A^, A , and A„ are shown in Fig. 3 at representative frequencies 
x y z 

common to the 14 cases. The averages of 14 values are plotted with error bars 
designating the one standard deviation error in the mean. (The point at 8. 3 MHz 
is an average of only seven cases. ) The amplification for the horizontal com- 
ponents starts near unity at the lowest frequency and rises significantly with 
frequency to values approaching 4 at the highest frequencies. For the vertical 
component the amplification remains near unity. 


The results in Fig. 3 were anticipated by theoretical analyses. Blank and 
Sill (1969) derived a model for the TE mode and Schubert and Schwartz (1969) 
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derived a model dealing consistently with the combined TE and TM modes in a 
wave field. In both treatments a thin current layer was assumed to surround the 
moon such that outside of the current layer the field had the free stream value, 
the normal component of which (since V • B = 0) was conserved across the layer. 
The approximation is appropriate on the sunlit side where the solar wind plpsma 

is moving into the moon at a speed both supersonic and superAlfvenic, so that 

/ 

in the magnetohydrodynamic approximation no information regarding the moon's 
presence may be conveyed upstream. The confinement permits very large ampli- 
fications of horizontal components. The current layer is not appropriate on the 
dark side; however, model calculations show that measurements on the sunlit 
side are not altered more than a few percent by the dark side portion of the 
assumed current system, and the symmetry is necessary for tractibility. 

Attenuation of the vertical component at a given frequency would be sub- 
stantially complete at a depth where the electrical skin depth for that frequency 
is small. If the moon were immersed in a vacuum, the attenuation would drop 

O _____ 

off as r and be observed at the surface. Tfte near unity values for the vertical 
amplification imply that a current layer is indeed present above the lunar surface 
and much closer than the underlying volume of high electrical conductivity. 

The scatter in the data has not yet been satisfactorily explained but may 
be due to plasma noise effects and departure from spherical symmetry. 
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The magnetic field amplification at the lunar surface, observed experi- 
mentally, is believed to give information about the electrical conductivity in 
the lunar interior. The moon is considered here to have an electrical con- 
ductivity o(r) varying only with radial distance. The magnetic field oscillation 
in the solar wind is defined as 



The cartesian coordinate system (£, rj, f) with unit vectors v_, £ is fixed 
relative to the moon, which moves through the solar wind in the negative ? 
direction with speed v relative to the wave front. (For transverse waves, 
neglecting propagation speeds, v is the aberrated solar wind speed.) The ampli 
tude is H Q , the wavelength X, and the frequency of the oscillation is f. In the 
lunar interior the solution of Maxwell’s equations for the TE mode is repre- 
sented by a potential SI satisfying 


V^S2 + k 2 S2 = 0 


where 


k 2 = 4n 2 f 2 ^e + i2nffur 
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and p and e are taken as free space values. The potential is found as a sum of 
product solutions of the form 


oo i 

£1 = sin <t> ^ F g (r)p g (cos 0) 

where r, 0, and 0 are spherical polar coordinates with f the polar axis and p* 
are associated Legendre polynomials. The boundary condition for the TE mode 
is the matching of the lunar surface normal component of the magnetic field to 
the free stream value, as discussed earlier, confining the lunar disturbance to 
the immediate proximity of the moon. 

Equations for the calculation of the fields are given elsewhere (Sonett et al. , 
197 Id); alternatively they may be written for the TE case by taking the solution 
to the earth induction problem (Lahiri and Price, 1939, or see Chapman and 
Bartels, 1962, Section 22. 13, equations 23-31) and making proper modifications 
for the changed boundary conditions. For the work reported here, the transfer 
function is defined as the ratio of the magnitudes of transverse components for 
the first mode of the modal expansion. 

Using the values of Ay and A z averaged from 14 swaths (of 1- and 2-hour 

duration) Sonett et al. have calculated a curve of amplification versus frequency. 

Because of the significant differences in Ay and A^, not yet explained, fits were 

— I 2 o il/2 

made separately for Ay, A z , and A = |0. 5 (Ay + A^m 
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The inversion method consists of starting with a conductivity profile, calcu- 
lating an A versus f curve and minimizing by the least squares method the de- 
viation of the calculated A with A, Ay, or A z . In order to do this the conductivity 

profile was characterized by eight parameters, namely the conductivity at 

r = 800, 1200, 1400, 1450, 1490, 1510, 1550, and 1740 km. The conductivity 
for 0 < r < 800 km was considered constant, the method not being very sensitive 
to conductivities at these depths. At all other values of r, log r is defined as 
a linear interpolation from neighboring log r values. In each iteration the ampli- 
fication for the model is calculated at eight frequencies, namely 0.83, 1.7, 5, 12, 
17, 22, 25, and 35 MHz, and the least squares fit at those frequencies between 
the model and the data is minimized by successive computer iterations using the 
Newton Raphson method. 

The method showed significant convergence in five iterations. Fig. 4 shows 
the resulting fit of a model amplification curve in which the rms value A was 
fitted. Fig. 5 shows the conductivity profile found by the routine that produces 
this result and profiles found by matching Ay or A z independently. In each case 
the profile is found to have a maximum near r = 1500 km. The conductivity starts 
at a low value at the surface as is expected because of the mean temperature of 
-30°C and the strong temperature dependence of candidate lunar materials. Be- 
cause the current system is effectively cut off, the method is imprecise at low 
conductivities. The conductivity then rises with depth to a maximum of about 
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6 x 10 ohm meter at r = 1500 km and then decreases by some 2 orders 
of magnitude in 100 km before resuming a rise with depth. 

The temperature profile can be inferred from the conductivity profile only 
indirectly, since it also depends on composition and conductivity-temperature 
functions. If candidate materials are ordered according to conductivity, one of 
the more conducting is a basalt such as the Lunar Sample 10024.22 analyzed by 
Nagata et al. (1970; see also Schwerer et al. , 1971). One of the less conducting 
is an olivine (see England et al. , 1968). The conductivity gradient from the 
surface to r = 1500 km is approximated by a 2°C/km temperature gradient for 
the basalt function and a 4°C/km gradient for the olivine, with most other candi- 
date materials lying somewhere in between. It is not likely that the temperature 
can decrease with depth between r = 1500 km and r = 1400 km because any such 
gradient occurring during the moon's formation would have smoothed out during 
the lunar history, and localized present-day heat sources at that depth are un- 
likely. If this reasoning is correct, the transition from 1500 to 1400 km must be 
accompanied by a modest increase in temperature as weir as a large decrease 
in conductivity. The transition is then plausible if the conductivity function at 
1500 km is relatively high, like that of a Nagata basalt, and the function at 
r = 1400 km is low, like that of an olivine. This reasoning leads to a moderate 
temperature profile such as sketched in Fig. 5, with a temperature on the order 
of 700°C at r = 800 km and approximately 800°C at the center. Since the 
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sensitivity of the model is low near the center, the central temperature is in- 
ferred from the fit of lunar thermal history models to the profile outside of the 
boundary r = 800 km (Fricker et al'. , 1967). The thermal profile curve is con- 
sistent with a lunar thermal history model with 25% chondritic radioactive con- 
centration; i.e. , one must assume that the moon had been relatively deficient in 
radioactive concentration compared to the chondritic meteorites during its for- 
mation to avoid higher temperatures in the present epoch, assuming traditional 
values for thermal conductivity. Acceptance of a solid state thermal conduction 
mechanism such as proposed by Runcorn (1962; see also Turcotte and Oxburgh, 
1969) would not alter the present-day thermal profile estimate but would allow 
the moon to have been hotter at an earlier time. 

In a parallel investigation using a different set of magnetic data Dyal and 
Parkin (1971) derived a monotonic lunar profile. The data were restricted to 
periods when the LSM was in the lunar night and consequently separated from 
the confining pressure of the incoming solar wind. Any lunar induced magnetic 
perturbation on the dark side is free to propagate into the essentially plasma- 
free cavity in the solar wind shadow. The geometry is complicated by the fact 
that the field direction is generally oblique to the axis of the shadow region. 

No unified model of the moon's reaction that includes both the sunward confine- 
ment and the cavity region has been developed. For traetibility, just as Sonett 
et al. (1971a, d) assumed a spherical current layer for the dayside analysis, 
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Dyal and Parkin assumed the moon to be immersed in a vacuum, since the dis- 
tance from the lunar measurement site to the confining currents is large. 


Dyal and Parkin analyzed the vertical component of many step function 
transients observed in the dark side. They fitted the data to a two-layer model 
moon of radius R^ with r = R-^ denoting the boundary between the inner layer 
of conductivity <r . and the outer nonconducting layer. Theory predicts that if 
the vacuum field surrounding the model moon is homogeneous and undergoes a 
unit step increase, the vertical component of the surface field rises as (Dyal 
and Parkin, 1971) 


B = 1 



jpf ex P 


/-S 2 TI 2 t 

VO "1 R ! 


The initial value at t < 0 is B = 0, and for t ,> 0, B = l-(Rj/R m )3. The final 

2 2 

value is unity. The dominant time constant is t = cr^ R^ /it , controlling ~60% 
of the series at the onset and dominating more and more as the higher order 
terms die out. 


In Fig. 6 is shown the data for one of the many steps analyzed by Dyal and 
Parkin; in this case the step was negative -going. The dashed line shows the fit 
of their model, with 0.95 R^ < R^, o- = 1. 5 x 10~ 4 ohm -1 meter - * and the time 
constant t ~ 55 seconds. The curve marked "step response" is the response 
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under the same conditions of a model similar to that of Fig. 5, a single current 
layer profile in which the conductivity is considered negligible except at r = 1505 
km where a layer of thickness Ar has a conductivity o characterized by oAr =70 
ohm . The initial fractional drop is 1-(R 1 /R m ) , for this model as well as for 
the two-layer model, and for R^ = 1505 km the initial drop does not conform to 
the data, as is shown in the figure. Schubert and Colburn (1971) suggest, how- 
ever, that the input function is imperfectly known but that since most transients 
in the solar wind are frozen into the plasma, the entire moon does not see the 
field change at the same time. A more appropriate concept is that of a wave 
front sweeping over the moon in a finite length of time. Schubert and Colburn 
approximate the effect of the moving wave front by a ramp function by which the 
surrounding field, assumed homogeneous, rises to its final value in 15 seconds. 
For this driving function the current layer model responds according to the curve 
labeled "ramp response", which is satisfactorily close to the data. 

The data of Fig. 6 represent but one case of the 10 used by Dyal and Parkin 
to obtain a size and conductivity for the core: 6 = 1.7 ±0.4 x 10~* ohm~ x meter _x , 
Rj > 0. 95 Rjjj. Consequently, while a fit to the data of Fig. 6 is not definitive 
Schubert and Colburn make two points: (1) it is important to consider the correct 
driving function for the moon, and (2) if the ramp approximation is valid, the 
profile of Fig. 5 and also two-layer profiles with R-^ < 0. 95 R m are eandiates for 
fitting the transient data. 
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Dyal and Parkin (1971) report a long tail on their transient response curves 

representing time decays longer than 4 minutes. The fit to these data implies a 

region of still higher conductivity deeper in the moon, represented in the model 

-2 -1 -1 

by a third layer of conductivity > 10 ohm meter and R 2 > 0. 6 
= 1040 km. The implications for the thermal profile are a higher inner tem- 
perature, as high as 1240°K = 967°C. 

It has been questioned whether the conductivity maximum in the conductivity 

profile is necessary to fit the data because a monotonic profile is simpler and 

would appear to follow the experimental data fairly closely. A comparison of 

the fits of various models is shown in Fig. 7. A two-layer model is shown, with 

-4 -1 -1 

a core of radius Rj = 1560 km and a conductivity a = 7.6 x 10 ohm meter . 

This model was attained as a best fit to the A data at the given radius and is close 

-4 -1 

to the model of Kuckes (1971) for which R 1 = 1580 km and a = 6 x 10 ohm 

meter *. The difference between the fits of the two-layer and the conductivity 

peak model is principally at the three lowest frequencies: a two-layer model 

fitting the middle range cannot dip sharply enough to approach the data at the 

lower frequencies. Other curves shown for comparison are those of model 3 

of Sill (1971) and Ness (1969). Other models investigated by Sill do not provide 

a better fit than his model 3. The model of Ness had R^ = 1426 km and 
-5 -1 -1 

o = 8 x 10 ohm meter . This model was based on an interpretation of a 
transient event observed by Explorer 35. Sonett, Mihalov and Ness (1971e) have 
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since proposed an alternate explanation for the data that would remove the basis 
for that model. It was concluded that Explorer 35 was too far from the moon to 
measure a lunar response. 

-2 

A difficulty with model fits for frequencies greater than 10 Hz lies in, the 
wavelength of the driving function (Schubert and Schwartz, 1971b). Convection 
in the solar wind dominates; wave propagation velocities in the frame of refer- 
ence of the solar wind plasma are generally much less than the bulk velocity of 
the solar wind so that the wavelength can be approximated by 

X = (n • v)/f 

where n is the unit normal to the wave front and v the solar wind bulk velocity. 

-2 

For frequencies below 10 Hz, x/Rjn » 1 for nearly all wave directions. Model 
fits by Sill and Kuckes are based on infinite wavelength. For the work reported 
here, n • v is taken to be 400 km sec 1 , introducing differences in the high fre- 
quency response and also in the definition of the driving function using a modal 
expansion. 

It is possible to extend the model to include the higher order spatial har- 
monics. The extended model, more accurate at high frequencies, has been pre- 
sented by Schubert and Schwartz (1971b). The use of the model at frequencies 
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_2 

above 10 Hz will depend on the feasibility of sorting out correct wave vector 
directions in the experimental data. 

Conclusion _ . _ . 

Measurement of a steady field at three places on the lunar surface implies 

j 

that the surface layer was magnetized during cooling by a field estimated to have 
been at least of 1000 gamma amplitude (Helsley, 1970). The source of this 
magnetizing field is not clear (cf. Sonett and Runcorn, 1971). Had a dynamo 
once been active in the lunar interior, using normal values of thermal con- 

/ 

ductivity, the moon would now possess a molten core. This appears unlikely 

from the electrical conductivity calculations, and also from rigidity arguments. 

\ 

If the field were an extension of the photospheric field, its value at the photo- 
sphere would be very large; under conditions where the field pressure dominated 

-3 5 

the solar wind, the field would decrease as r , requiring 10 gauss at the 

photosphere. If the solar wind had sufficient momentum density to shape the 

-2 

magnetic field, the radial field would follow an r law (Parker, 1963). The 
field requirement would then be more modest, 600 gauss, but the solar wind 
density would be orders of magnitude greater than the present day value, a 
situation not expected to have occurred after a possible T-Tauri stage of solar 
evolution (Sonett et al. , 1968). Polarity reversals due to the sun's rotation 
would also make it unlikely that the sun's field could magnetize a rock layer 
slowly cooling through the Curie point. 
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Immersion of the moon in the earth's field to attain a 1000 gamma magnetizing 
field implies a synchronously rotating moon at an orbital distance dangerously 
close to the Roche limit, unless the earth's field had been far greater than present 

day values (Dyal et al. , 1970a). _ A satisfactory explanation for the required 

magnetizing field has not yet been found. 

The conductivity profiles obtained by the dayside data imply a discontinuity 

at r = 1500 km. The conductivity maximum implies a change at that depth from 

a more to a less electrically conducting material. Alternatively the presence at 

that depth of a relatively thin highly conducting layer as suggested by Urey (see 

Urey et al. , 1971) would also fit the electrical conductivity data. In this case the 

2 -1 

oAr product would be -10 ohm and the lunar interior temperature might be 
lower than suggested if a rock other than olivine is postulated for the deep interior . 

The barrier at a depth of some 250 km appears even for models with mono- 
tonic conductivity profiles, although there is some question from the darkside 
transient data as to whether a shallower depth is more appropriate. A barrier 
implies a change of composition or state at that depth, possibly the lower limit 
of a surface melting process. The process may be tied to the surface melting 
required to explain the ages of lunar samples. Reynolds et al. , (1971) have 
considered models for the surface heating process by a combination of accretional 
and radioactive means. They find certain combinations that fit the timetable of 
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lunar rock. ages. Wright (1971) suggests that the conductivity profile can be 
explained by oxygen depletion. 

Except for the Explorer 35 survey the positive magnetic measurements of 
the moon have been confined to three locations for the steady field and one loca- 
tion for the time varying field. The successful deployment of the Apollo 15 LSM 
in August 1971 promises additional detail, and together with the Apollo 16 LSM 
will allow consideration of possible lunar asymmetry (Schubert and Schwartz, 
1971a). Other refinements to be investigated are the separation of the TM and 
TE modes and the determination of the lunar response for various k vectors. 
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FIGURE CAPTIONS 

Fig. 1. Steady magnetic field vector at three locations on the moon. Site 12 
is the Apollo ALSEP location; sites 14A and C’ are the two measuring sites on 
the Apollo 14 moonwalk. The maximum angle between any two vectors is 84°. 

(From Dyal et al^, 1971) _ 

Fig. 2. Power spectral density of the interplanetary field near the moon 
(Explorer 35) and the lunar surface field (Apollo 12) for the same 2-hour time 
period. The two horizontal components are shown: Z is north and Y is east at 

the Apollo 12 site. Other 2-hour segments show variations from this example 
but all show an amplification in the horizontal components increasing with frequency. 

Fig. 3. Amplification, the ratio of amplitudes of lunar surface and free stream 
magnetic fields, shown as a function of frequency. The amplification is defined 
as the square root of the ratio of power spectral densities. The error bars for 
frequency are the windows defined by the lags in the autocorrelation calculation. 

The error bars in amplification are the one standard deviation limits determined 
from the means of 14 data spectra. Amplification occurs for the north and east 
components but remains near unity for the vertical component. 

Fig. 4. Amplification as a function of frequency. The data are as in Fig. 3, with 
the tangential amplifications combined by the relation A = [o. 5(A^ + A^)] 1//2 . The 
solid line is the amplification calculated from a model described in the text. The 
amplification of the model is fitted to A values at frequencies of 0. 83, 1.7, 5, 12, 


17, 22, 25, and 35 MHz. 
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Fig. 5. Lunar electrical conductivity as a function of radial distance, according 
to calculations based on the amplifications Ay, A^, and A. As depth increases 
the conductivity rises to a peak at r = 1500 km, decreases to a minimum at 
r = 1400 km and then continues to rise. Six other cases were run, using the one 
standard deviation limits of the As for calculating the conductivities; these all 
contain the peak and lie in the shaded area shown. Plotted also is a tentative 
version of a lunar thermal profile. The assumption of a relatively smooth 
thermal profile and known conductivity functions requires a core conductivity 
similar to that of olivine with a transition to a more highly conducting material 
at r = 1500 km. 

Fig. 6. Transient observed in the vertical component of lunar surface magne- 
tometer data while the magnetometer site was in lunar night. Eddy currents in 
the interior inhibit rapid changes in this component. The dashed line is the re- 
sponse of the Dyal and Parkin model obtained from many of these transients 
(Dyal and Parkin, 1971). Another model, described in the text, has the re- 
sponses shown depending on whether the input function is best described as a 
step or a ramp (Schubert and Colburn, 1971). 

Fig. 7. Amplification as a function of frequency for several lunar electrical 
conductivity models. The data and the curve labeled "best fit" are as shown in 
Fig. 4. The curve labeled "two-layer" is a best fit under the constraints of a 
two-layer model when = 1560 km. The conductivity peak in the model is 
necessary in order to come close enough to the three lowest frequency data 


points. 
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